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1 .  Introduction 

In  this  note  we  present  improvements  in  the  methods  for  treating  the 
Compton  current  and  the  air  conductivity,  for  use  in  EMP  calculations. 
Specifically,  we  shall  give: 

(a)  a  method  of  treating  the  effoct  of  multiple  scattering,  as  well 
as  slowing  down,  of  the  Compton  recoil  electrons; 

(b)  a  method  of  treating  the  time  lag  between  the  production  of  the 
primary  ionization  and  of  the  total  ionization; 

(c)  a  discussion  of  the  problem  of  equilibration  of  the  free  elec¬ 
trons,  from  the  initial  ionization  spectrum  to  the  spectrum 
appropriate  to  the  instantaneous  value  of  E/p  . 

2.  Slowing  Down  of  Compton  Recoil  Electrons 

The  energy  loss  of  relativistic  electrons  in  moving  through  matter  is 
given  by  Bethe's  formula  (Ref.  I);  per  unit  track  length,  it  is 

-  «  M2|l  [i]  (1) 

ds  mv2  1  J  » 
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In  these  formulae,  N  is  the  density  of  atoms,  Z  the  atomic  number, 
e  ,  m  and  v  the  electron  charge,  rest  mass,  and  velocity,  c  is  the 
velocity  of  light,  and 


0 


v 
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/l-82 

For  air  we  take  the  effective  atomic  number 

Z  =  7.2 


(3) 


(4) 


JS 

% 
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and  the  mean  excitation  potential  I  (Ref.  1}  to  be 


I 


80.5  ev 


(5) 


l»?  solving  the  equations  of  motion  of  a  Compton  recoil  electron,  one 
can  take  account  of  the  energy  loss  by  imagining  that  the  electron  experi- 
eiices  a  force  directed  opposite  to  v  ,  of  magnitude  given  by  the  right- 
hand  side  of  £q.  (1).  This  procedure,  however,  does  not  yet  take  into 
account  the  multiple  scattering  of  the  Compton  electrons,  which  has  the 
following  effect.  A  beam  of  electrons,  directed  initially  in  the  x-direction, 
is  greatly  broadened  in  angle  by  the  scattering,  so  that  the  electrons  spend 
a  lot  of  their  energy  moving  in  directions  oblique  to  the  x-direction.  Eq. 

(1)  is  correct  for  the  energy  loss  per  element  of  length  of  the  actual  trajec¬ 
tory  of  an  individual  electron,  but  its  use  with  ds  replaced  by  dx  will 
overestimate  the  average  distance  electrons  will  travel  in  the  x-direction 
before  stopping. 

rn  the  p3st ,  we  have  taken  account  of  the  multiple  scattering  by  not¬ 
ing  that,  experimentally,  the  mean  range  of  electrons  in  the  Mev  region  is 
about  2/3  of  the  total  (or  "extreme'')  range.  We  have  therefore  simply 
multiplied  the  drag  force  (right-hand  side  of  Eq.  (1))  by  a  factor  3/2. 

This  procedure,  however,  gives  a  wrong  distribution  of  ionization  along  the 
path  of  the  Compton  electron.  The  ionization  per  unit  length  is  proportional 


f 


l‘i  it ■.**» 


to  "  »  am*  is  important  because  it  determines  the  air  conductivity. 

Now  at  the  beginning  of  a  Compton  electron  track,  the  multiple  scattering 

dW 

has  not  yet  had  much  effect,  and  multiplying  by  3/2  overestimates 

the  ionization  per  unit  length.  Similarly,  the  ionization  density  is 
underestimated  near  the  end  of  the  range.  In  EMP  from  high  altitude  bursts, 
the  peak  electric  field  tends  to  come  when  the  Compton  electrons  are  near 
the  beginning  of  their  range;  an  overestimate  of  air  conductivity  by  a 
factor  3/2  may  lead  to  an  underestimate  of  the  peak  electric  field  by  the 
same  factor. 


Obviously,  we  need  a  treatment  in  which  the  multiple  scattering  is 
allowed  to  build  up  along  the  trajectory. 

3.  Multiple  Scattering  of  Compton  Electrons 

We  shall  use  the  simple  fora  of  multiple  scattering  theory  due  to  Williams 
(Ref.  2),  which  will  be  adequate  for  our  purpose.  According  to  this  theory, 
the  mean  squared  angle  W  of  an  initially  collimated  beam  increases  with 
distance  traveled  as 


dW  _  SttNzV  ryl 

ds  ~  y2m2vi'  1  ^ 


(6) 


where 

[2]  =  £n  (f~)  .  (7) 

\  min/ 

Here  is  the  minimum  angle  of  scattering,  below  which  the  scattered 

angular  distribution  falls  substantially  below  that  given  by  the  Rutherford 
formulae.  For  scattering  by  neutral  atoms,  as  in  our  case,  8  is  deter¬ 
mined  by  the  screening  of  the  nuclear  charge  by  the  atomic  electrons.  For 
the  Fermi-Thom&s  atomic  model,  Mott  and  Massey  (Ref.  2)  give 
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In  Eq’s.  (6)  and  (7)  ve  have  departed  slightly  from  standard  theory, 
in  which  02  is  the  mean  square  angle  in  only  the  gaussian  part  of  the 
angular  distribution.  Then  the  argument  of  the  logarithm  in  Eq.  (7) 

.comes  (£ij/®min)  ,  where  0j  depends  on  s  ,  and  is  the  angle  below 
which  the  scattering  is  multiple.  However,  we  are  using  82  as  the  mean 
square  angle  of  the  entire  angular  distribution,  in  which  case  the  argument 
(2/6n-n)  is  correct  (Ref.  3). 

Eq.  (6)  is  strictly  correct  only  when  9^  is  small,  where  successive 
scattering  angles  add  linearly.  Consider,  however,  the  obliquity  factor 


1 

cos  8 


This  factor  relates  the  differential  track  length  ds  to  the  distance  dx 
in  the  direction  of  initial  collimation. 


For  small  9  , 


so  that  Eq.  (6)  becomes 


4T?rzV 

yTmspr~ 


For  large  mean  angles  of  scattering,  this  equation  has  the  virtue  that  near 


A 
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the  end  of  the  range,  where  n  will  become  infinite,  the  mean  cos  8 
will  become  zero,  according  to  Eq.  (10);  thus  all  correlation  to  the 
initial  direction  is  lost.  This  is  a  desirable  result,  whereas  having 
Q7  become  infinite  (as  results  from  Eq.  (6))  is  not.  Ke  shall  there¬ 
fore  use  Eq.  (13)  instead  of  Eq.  (6),  and  note  that  n  starts  at  unity 
for  an  initially  collimated  beam. 

4.  Mean  Range  cf  a  Collimated  Beam 

Suppose  that  we  have  a  monoenergetic  beam,  collimated  in  the  x 
direction,  entering  a  slab  of  material.  Then  the  equations  for  the  change 
in  kinetic  enegv  K  and  obliquity  factor  n  (we  drop  the  bar  on  n)  are 


an 

dx 

f(fc’)  n 

(14) 

dW 

dX 

•g(h')  n 

(15) 

where 

f(W) 

4irNZ2eV  r-ii 

(mc*)2^2-!)2  1  J 

(16) 

gOO 

-  ^rNZeV  [13  . 

mcz(y2-l)  1  J 

(17) 

Dividing  Eq. 

(14)  by  Eq.  (IS)  and  using 

W 

(Y-J)  me2 

(18) 

we  find 

an  = 

.  -11—  ill 
y2-i  m  { 

(19) 

Now,  it  turns  out  that  the  ratio 


r 


It} 


(20) 


is  very  nearly  independent  of  energy  W  or  y  .  For  example,  for  air  (see 
Eq's  (4)  and  (5))  we  find: 


K 

y 

r 

51  Kev 

1.1 

0 . 262 

255  Kev 

1.5 

0.266 

511  Kev 

2.0 

0.269 

1.02  Mev 

3.0 

0.270 

1.53  Mev 

4.0 

0.271 

Clearly,  for  our  purposes,  T  can  be  taken  as  constant.  He  shall  use  the 
values. 


F  =  0.269  in  air  , 

T  =  0.276  in  aluminum 

'T  =  ISOev  for  aluminum) 


(21) 


Kith  the  approximation  of  constant  F  ,  Eq.  (19)  can  be  integrated 
directly,  with  the  result. 


r\ 


1  +  zr  In 


(Y o - 1 )  (>■*•!) 
l(Yo*l)(Y-l) 


(22) 


Here  y0  is  the  value  of  y  for  the  initial  energy  of  the  electrons  in 
the  beam. 


On  using  Eq.  (22)  in  Eq.  (15),  one  can  integrate  (numerically)  to 
find  the  mean  range  including  the  effect  of  multiple  scattering.  We  have 
done  this  for  aluminum,  and  also  have  calculated  the  extreme  range,  which 
is  obtained  by  omitting  the  factor  n  in  Eq.  (15).  In  Fig.  (1)  we  compare 
our  results  with  data  of  Marshall  anri  Ward  (Ref.  4)  on  the  transmission  of 
electrons  through  foils  of  aluminum.  In  the  figure,  the  solid  curves 
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represent  Marshall's  and  Ward's  measurements  rf  the  fraction  of  initial 
electrons  transmitted  as  a  function  of  thickness  of  the  f^il,  for  various 
starting  electron  energies  from  0.32  to  1.69  Mev.  The  sqc'fts  placed  on 
the  abscissa  are  at  the  extreme  ranges  computed  by  vs  usir;-  the 
formula  directly;  the  agreement  with  the  end  points  of  the  e-xperiscntal 
curves  is  ver>'  good.  'Ti'*  circles  placed  on  the  experimental  curves  are  at 
the  mean  ranges  computed  by  .j  using  F.q's  (IS)  and  (22).  Our  theory  would 
be  well  verified  if  the  circles  also  represented  the  mean  abscissae  of  the 
experimental  curves.  We  consider  the  agreement  between  theory'  and  data  to 
be  quite  adequate  for  our  purposes. 

In  Sec.  5  we  shall  combine  this  treatment  of  energy  loss  and  multiple 
scattering  with  equations  of  motion  for  the  electrons  in  the  presence  of 
fields.  Before  turning  to  that  task,  we  record  in  Table  I  the  mean  and 
extreme  ranges  of  electrons  in  air,  as  computed  from  the  above  theory, 
since  these  results  arc  of  general  interest. 

Tabic  i 

Ranges  of  Electrons  in  Air 


Kinetic  Energy 

W,  Mev 

Mean  Range 

Rfij>  gm/cra2 

Fit  to 

Kean  Range 

Extreme  Range 

Rp,  gm/cm2 

s 

1 

j 

* 

0.05 

0.003G 

0.0029 

0.0045 

:j 

if 

0.1 

0.0036 

0.0100 

0.0157 

£ 

i 

0.2 

0.030 

0.032 

0.050 

-9 

0.5 

0.120 

0.125 

0.198 

T 

-1 

1.0 

0.304 

0.303 

0.489 

i 

2.0 

0.711 

0.695 

1.085 

J 

5.0 

2.01 

1.89 

2.78 

£7 

It  may  be  useful  to  note  that  the  mean  range  can  be  fitted,  over  the 
energy  range  in  the  table,  by  the  simple  formula 

„  ,  ,  2,  0.40  H2 

RflCgm/c®  )  =  06^?"  M 

Table  I  also  contains  values  computed  from  this  formula. 
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5 .  Equations  of  Motion  in  Presence  of  Electromagnetic  Fields 

We  shall  now  write  equations  of  motion  for  the  Compton  recoil  electrons, 

-y- 

in  the  presence  of  electric  and  magnetic  fields  E  and  B  .  and  taking 
account  of  energy  loss  and  multiple  scattering. 

If  no  account  were  taken  of  multiple  scattering,  but  energy  loss  were 
included,  the  equations  of  motion  would  be 


M  I?  +  Jr*  s)  -  sf>«)  lie 


Here  g(W)  is  the  energy  loss  function  defined  by  Eq.  (17),  r  is  the 

position  of  the  electron,  t  is  the  time,  and  p  is  the  electron  momentum. 

")■ 

The  relations  between  p  and  v  and  W  are 


.£ _ 


/pJx4+m4c"£ 


Integration  of  Eq's  (24)  and  (25)  would  lead  to  precise  tra jectories,  with 
the  electrons  being  deflected  by  the  fields  and  slowed  by  the  energy  loss 


term. 


Now  consider  the  multiple  scattering.  Since  wo  have  already  accounted 
for  energy  loss,  it  is  necessary  only  to  account  for  the  fluctuations  in 
the  angles  of  p  ,  without  additional  changes  in  the  magnitude  of  p  . 

“v* 

The  result  of  the  many  small  and  random  changes  6p  (with  5p  1  p)  will 
be  to  make  a  probability  distribution  in  p  .  The  value  of  p  calculated 

*k 

from  Eq.  (24),  without  scattering,  will  be  the  central  p  of  this  distribu¬ 
tion,  and  the  distribution  will  tend  to  be  over  momenta  having  the  same 
magnitude,  as  in  Fig.  (2a),  although  the  electric  field  ~£  car.  le.td  to 
situations  as  illustrated  in  Fig.  (2b). 


1  -'TrCa'iT  aa.>_  * 
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(a)  (b) 

Figure  2.  Relation  of  central  momentum  p 
to  distribution  at  a  given  time. 

Therefore,  we  retain  Eq's  (24)  and  (26)  for  the  central  momentum. 

The  average  velocity  of  the  distribution,  however,  is  less  than  the  value 

~r 

appropriate  to  p  ,  by  the  mean  obliquity  factor  n  .  Therefore,  we 
replace  Eq.  (25)  by 


,**•  -► 

dr  v 

dt  *  n  (27) 

-b 

where  v  is  calculate'  from  Eq.  (26)  from  the  central  momentum.  The  equa¬ 
tion  for  n  is  taken  over  from  Eq.  (14), 


dn 

dt 


v  f(W) 


(28) 


Repeating,  we  take  the  equations  of  motion  to  be  Eq's  (24),  (26),  (27) 
and  (28).  The  Comptcn  current  contributed  by  an  electron  (really  a  distribu 
tion  of  electrons)  is  taken  as 


_e__  dr  __  J_eJ  v 

c  dt  ~  c  n 


(29) 


The  time  rate  of  production  of  ion  pairs  by  the  electron  is  taken  as 


|i  =  v  g (W)/Wj  .  (30) 

where  Wj  is  the  energy  expended  per  i" -  pair. 

Notice  that  the  equations  of  motion  chosen  above  reduce  to  Eq's  (14) 
and  (15)  in  the  absence  of  fields.  For  in  that  case  the  central  momentum 
will  remain  in  the  x  direction  if  initially  pointed  in  it,  and  from  Eq. 
(27)  we  find. 


Thus  Eq.  (28)  immediately  reduces  to  Eq.  (14).  Further,  since 


r  -r 

v  •  dp 


Eq.  (24)  reduces  to  Eq.  (15). 

In  some  cases,  it  is  convenient  to  use  the  retarded  time  (r  =  distance 
from  burst  point). 


t  -  - 
c 


rather  than  t  iri  the  equations  of  motion.  The  change  to  retarded  time  is 
accomplished  through  the  relation 


(i  -  y~) « 

\  nc  / 


Here  v^  is  the  radial  component  of  the  velocity  computed  from  Eq.  (26). 
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6. 


Change  in  n  Due  to  Fields 


Eq.  (28)  shows  how  n  increases  with  time  due  to  the  multiple  scatter 
ing.  We  shall  see  now  that  additional  changes  in  r)  occur  if  the  electric 
field  has  a  component  E„  parallel  to  the  central  p  . 

First,  it  is  clear  that  the  magnetic  field,  and  the  component  Ej_  of 
the  electric  field  perpendicular  to  J  cause  no  direct  change  in  n  .  In 
Figure  3  let  the  solid  vectors  represent  central  and  extreme  momenta 


Figure  (3)  Change  in  momenta  induced  by  a  force,  (a)  perpendi¬ 
cular  to  p,  and  (b)  antiporaiiel  to  p. 


initially,  and  the  dashed  vectors  represent  the  same  after  action  of  a 
perpendicular  force  in  (a)  and  an  antiparallel  force  in  (b).  In  case  (a) 
the  entire  bundle  of  momenta  is  deflected,  but  the  argular  spread  is  not 
changed  in  first  order.  But  in  case  (b)  the  angular  spread  obviously 
increases.  If  the  force  were  parallel  to  p  ,  the  angular  spread  would 
decrease. 

For  a  given  momentum  in  the  bundle  of  momenta,  let  pj|  be  the  com¬ 
ponent  parallel  to  the  central  momentum,  and  p±  the  perpendicular  com¬ 
ponent.  The  angle  0  this  momentum  makes  with  the  central  momentum  is 
given  by 

tan  0  =  P|/P„  •  (35) 
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From  Oq.  (10)  one  can  find,  by  trigonometry. 


tan  9 


✓rjM 


Therefore, 


W 


Let  there  be  an  electric  field  E|(  parallel  to  the  central  p  . 
equations  of  motion  are 


Then  the 


-iel  H, 


lhen,  differentiating  Eq.  (37)  with  respect  to  time  yields 


Pl>,3,  ^ 
n  dt 


n2-i 
n  Pn 


iel  Cn  -l) 


Z-i'l  £  ' 


where,  in  the  last  line,  p  is  the  central  p  . 

Therefore,  to  include  the  effects  of  both  multiple  scattering  and  E„  , 
Eq.  (28)  must  be  replaced  by 


,2  ' 


v  f (W)  +  i e {  (n  -l)£— J 


7 .  Spatial  Spreading  of  a  Beam 

As  a  result  of  the  angular  spread  in  velocity  induced  by  multiple 
scattering,  an  initially  collimated  beam  will  also  spread  in  space,  around 
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via 


the  central  point  determined  from  the  equations  of  motion  derived  above. 
We  shall  discuss  this  spreading  in  an  approximate  way,  valid  when  the 
angular  spread  is  small,  and  in  the  absence  of  fields.  We  shall  use  the 
Boltzmann  equation,  which  determines  the  evolution  of  the  distribution 

•V  — 

function  $(r,v)  in  phase  space, 


i*  ♦  v  .  v« 

Q  t 


collision  terms 


The  collision  terms  have  two  effects.  First,  the  magnitude  of  the 
velocity  v  of  a  given  particle  is  gradually  decreased  with  time.  The 
average  decrease  in  v  corresponds  to  the  energy  loss  given  by  Eq.  (1). 

There  are  also  fluctuations  in  the  energy  loss,  which  give  rise  to  a  spread 
in  the  magnitude  of  v  at  a  given  time;  this  in  turn  gives  rise  to  a 
spread  in  position  of  the  particle  at  a  given  time.  However,  this  spread 
in  position  is  less  important  than  that  induced  by  fluctuations  in  the 
direction  of  v  ,  which  we  discuss  below.  Therefore  we  shall  assume  that 
the  magnitude  v  decreases  smoothly  with  time,  and  ignore  the  spread  ir 
v  .  However,  we  shall  give  below  an  estimate  of  the  effect  of  the  spreau 
in  v  . 

The  other  effect  of  the  collision  terms  is  to  cause  fluctuations  in  the 
direction  of  v  .  Because  the  scattering  angle  per  collision  is  predominantly 

•V 

very  small,  the  spreading  in  direction  of  v  can  be  treated  in  the  diffusion 
approximation,  with 


collision  terms  -  a  7^  $  .  (42) 

Here  VQ2  is  ?.he  angular  part  of  the  Laplacian  operator  in  velocity  (rather 
than  coordinate)  space,  and  «  is  (one  fourth  of)  the  mean  square  scattering 
angle  per  unit  time.  This  quantity  has  already  been  used  in  Sec.'s  (3)  and 
(4)  above,  and  in  fact  we  have  (as  will  be  verified  below) 


a  (sec-1) 


7fW 


We  now  choose  n  Cartesian  coordinate  system  in  which  the  x-axis  is 
parallel  to  the  initial  velocity  of  the  collimated  beam.  We  introduce 
two-dimensional  vectors  for  the  two  othogonal  directions,  defining 


* 

r 


(y,z) 


fi 

4 


(44  J 


For  the  velocity,  8  will  denote  the  angle  between  v  and  the  x-axis. 

When  0  is  small,  it  can  be  thought  of  as  being  composed  of  a  deflection 
6  in  the  y-direction  and  a  deflection  6  in  the  2-direction.  (Two  num- 

y  z 

bets  are  required  to  specify  the  direction  of  the  velocity.)  Thus  we  may 
regard  6  as  a  vector. 


5  =  cey,y 


rr  2 
v  0 


32  32 

30^*  397 


(45) 


The  velocity  then  has  components,  to  first  order 
vx  =  v  (l  -  Y~j  parallel  to  x-axis 

v  =  v  ?  parallel  to  yz  plane 


(46) 


(Note  that  v  here  is  two  dimensional,  and  will  be  until  further  notice.) 


-J 


In  this  notation,  the  iioltzmann  equation  becomes 

^^(x^.Q.t)  ♦  v  (l  -  |1)  ||+  v$  •  V*  =  a702f  (47) 

In  this  equation,  v  is  understood  to  decrease  smoothly  with  time;  a  also 
changes  with  time,  from  Eq.  (43).  W’e  shall  normalize  the  initial  -p  to 
unity, 

yrp(x,v,t,0)dT  =1  ,  dx  =  dxdydzd0yd&2  -  (48) 

By  integrating  lq.  (47)  over  phase  space,  it  follows  that  ip  is  so  normal¬ 
ized  at  all  times, 
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We  have  assumed  here  that  x p  and  its  derivatives  vanish  fcr  large  values 
of  its  arguments. 

The  average  value  of  x  at  time  t  is 


f  x  tp(t)dT 


Multiplying  Eq.  (47)  by  x  and  integrating  over  phase  space,  we  find 


A  - 
dt  X 


*/(■  -  P) 


(>-r) 


Similarly,  multiplying  Eq,  (47)  by  62  and  integrating,  we  find 


at 


4o  (V0262  =  4) 


Comparing  this  result  with  Eq.  (28),  we  see  that  Eq.  (43)  is  verified. 


Let  us  now  calculate  the  spread  in  x  .  Kote  that 


(x-x): 


Thus  we  need  to  calculate  P  .  Multiplying  Eq.  (47)  by  x2  and  integrating, 
we  find 


x2 

dt 


0  ~  x92 
2v  x — 


Thus  we  need  to  find  x62  .  In  the  usual  way,  we  find 


c7  x0i 


v  9‘^i  -  +  4ax 


7T  P  .  — 
v  6  -  Tj—  +  4  ax 
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Now  we  need  6*  ,  which  is  determined  by 

~  0*  16  aF 

at 

These  relations  may  be  condensed  a  hit  by  noting  that 

4r  Ax2  =  x3*  -  2x  —•  =  v[x  82  -  x02  ] 


Tlie  quantities  on  the  right  here  have  to  be  found  from  Eq's  (.'>') .  (52). 
(55),  and  (56).  All  of  these  equations  have  to  be  Integrated  taking  into 
account  the  dependence  of  v  and  a  on  t  .  It  is  instructive,  however, 
to  obtain  an  approximate  answer  by  regarding  v  and  a  as  constants. 

One  then  finds 


e2 

4at  , 

X  = 

vt(l-at) 

6“  = 

52 (at) 2 

xe7  = 

vtjdat  - 

Ax2  = 

|  (vt)2 

4  (vt  -  x)‘ 

O 


From  Eq.  (62)  we  see  that  the  spread  in  x  is  conpanole  with  the  amount 
by  which  x  falls  behind  vt  ,  the  position  the  particle  would  have  if 
there  were  no  scattering.  This  result  is  ip  substantial  agreement  with 
the  curves  of  Fig.  (1). 

We  are  also  interested  in  the  lateral  spread  of  the  beam.  To  this 
end,  one  calculates  from  Eq.  (47), 

~  r5'  =  2  v  $*r  ,  (65) 

—  ?*r  =  vF  .  (64) 


In  the  approximation  that  v  and  a  arc  constants,  these  equations 
integrate  to 


$*r  =  2(vt) (at)  ,  (6S) 

rJ  =  2(vt)2(at)  ,  (66) 

F  =  U  =  (vt)2(at)  .  (67) 

One  sees  that  tne  lateral  spread  builds  up  more  rapidly  than  the  spread  in 
the  x-direction. 

If  we  had  considered  the  fluctuations  in  the  magnitude  of  v  ,  there 
would  have  been  an  additional  term  in 

=  (vt)2  (St)  ,  (68) 

the  mean  square  (fractional)  fluctuation  in  v  per  unit  time 

=  2  v2  St  (69) 

is  smaller  than  the  rate  a  by  the  ratio 


©(fix7) 

where  8  is 


The  rate  6 


8  a  _ l_ 

a  "  y2Z[2J 


(70) 


Here  [2]  is  given  by  Eq.  (9).  The  factor  y2  enters  because,  when 
v  -  c  ,  substantial  changes  in  energy  lead  to  only  small  changes  in  the 
velocity.  The  factor  2  occurs  because  energy  loss  is  due  to  collisions 
with  electrons,  while  angular  scattering  is  due  mainly  to  collisions  with 
the  nuclei.  The  factor  (2]  occurs  because  the  average  of  82  involves 
logorithmic  integral  in  Rutherford  scattering  while  the  average  of 
(6v}?  'v  0U  does  not.  In  air,  8  is  only  a  few  percent  of  a  .  Thus 


SO 


. 


fc- 


I 


while  the  contribution  to  Ax5”  in  Eq.  (68)  rises  as  a  lower  power  of  time 
than  that  from  Eq.  (62),  by  the  time  at  which  the  multiple  scattering 
becomes  important,  the  terms  in  Eq.  (62)  are  dominant. 


As  a  numerical  example,  consider  an  electron  beam  that  starts  with 
1  Mev  (ye  ~  3)  ,  and  let  us  examine  the  distribution  when  the  energy  has 
fallen  to  0.5  Mev  (y  =  2)  .  According  i-j  Eq.  (22),  we  will  then  have 
H  -  1.79  ,  for  air.  Then  from  Eq.  (10),  6^  =  1  ,  so  that  from  Eq.  (53), 
at  ~  1/4  .  From  Eq's  (39),  (62)  and  (67)  we  then  find 


Thus  the  lateral  spread  in  position  of  the  particle  is  about  2/5  of  the 
mean  distance  the  particles  have  traveled  at  this  time. 


We  see  that  the  spread  in  pc  /»  is  appreciable,  and  could  have 

noticeable  effects  if  included  in  calculations.  Ke  do  not.  propose 

that  the  spatial  spread  be  included  in  all  JMP  calculations,  but  believe 
that  the  method  of  Sections  (5)  and  (6),  using  the  average  position,  is 
adequate  for  most  purposes. 

However,  occasionally  we  may  want  to  make  a  calculation  including  the 
full  effects  of  multiple  scattering.  Ke  believe  that  the  best  way  to  do 
this  is  to  add  random  velocity  changes  6v  as  the  equations  of  motion  of 
the  Compton  electrons  are  being  solved.  Then,  of  course,  the  factors  n 
should  be  removed  from  the  equations  of  motion.  The  drag  force,  which 
gives  the  mean  energy  loss,  should  probably  be  left  in  the  equations  of 
motion.  Then  to  first  order,  the  fluctuations  Sv  should  be  such  as  to 
leave  the  magnitude  v  unaltered;  however,  fluctuations  in  6v  could 
also  be  included,  provided  the  average  6v  is  kept  equal  to  zero. 
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It  is  not  necessary  to  make  the  probrHi'  v  distribution  of  ov 
match  the  Rutherford  scattering  law.  As  1  s  jav[  «  v  ,  it  is 

necessary'  only  to  make  60^  and  ov2”  have  the  correct  mean  rates  of 
increase  (see,  e.g.,  Eq‘s  (52)  and  43))- 


8.  The  Ionization 

We  now  turn  to  another  effect  of  the  Compton  electrons,  namely  the 
ionization  they  produce  in  the  air.  This  is  important  because  it  deter¬ 
mines  the  electrical  conductivity  of  the  air. 


m. 


It  is  well  known  that  fast  electrons,  in  moving  through  air,  produce 
one  ion  pair  for  (about)  each  34  ev  of  energy  lest  by  the  fast  electrons. 

Not  all  of  this  ionization,  however,  is  produced  directly  by  the  fast 
electron.  Rather,  some  of  the  electrons  dislodged  by  the  fast  electron 
have  enough  energy  to  produce  further  ionization.  Ionization  produced 
directly  by  the  fast  electron  is  called  the  primary  ionization,  and 
further  ionization  is  called  the  secondary  ionization.  The  secondary' 
ionization  builds  up  gradually  in  time  after  production  of  the  primary 
ionization.  The  time  lag  in  formation  of  the  secondary'  ionization  is 
important  for  EMP  at  high  altitudes. 

The  cross  section  for  ionization  of  atoms  by  fast  electrons  was  also 
computed  by  Betbe  (Reference  5) .  The  cross  section  has  also  been  measured 
for  some  gases,  including  N2  and  02  ,  by  Schrara  and  collaborators 
(Reference  6).  Using  the  experimental  data  to  evaluate  uncertain  parameters 
in  Bethe's  formula,  we  have  deduced  the  following  formula  for  the  ioniza¬ 
tion  cross  section  per  atom  of  air: 


a.  (air  atom)  =  M,2  [3} 

1  RV  ' 


(71) 


where  is  the  Bohr  radius  and 


U]  .  a  a. 


+  2fny  -  1  *  —i 
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(72) 


with 


M  *  =  4.05,  V  ::  16  ev  for  air. 


(73) 


The  number  of  ion  pairs  I  produced  per  unit  length  by  the  fast  electron  is 


dI  -  v,, 

ds 


(74) 


where  N  is  the  density  of  air  atoms,  as  befcrrn.  Comparing  this  equation 
with  Eq.  (1),  we  can  compute  the  energy  lost  per  ion  pair  of  primary  ioniza- 
t  ion, 


dW 

dl 


SJ  2 

0 


Ml 

(3) 


(75) 


h'e  have  evaluated  this  result  for  several  electron  energies,  and  obtain  the 
numbers  in  Table  2. 


Y 

ev 


ion  pair 


1.1 

80.0 


Table  2 


1.5 

84.5 


2.0 

85.4 


3.0 

86.4 


4.0 

86.9 


It  is  seen  that  4y-  does  not  change  much  over  the  range  of  interest;  we 
shall  use  the  constant  value  86  ev/ion  pair. 


Thus  we  find  that  the  total  ionization  is  a  factor  86/54  =  2.52  larger 
than  the  primary  ionization,  i.e.,  each  secondary  electron  must  produce, 
directly  and  indirectly,  1.52  additional  electrons,  on  the  average.  W c 
need  to  find  out  how  this  residual  ionization  builds  up  in  tine. 

First,  it  may  occur  to  the  reader  that  the  Auger  effect,  following  ejec¬ 
tion  of  a  K-shell  electron,  should  result  in  a  practically  instantaneous 
increase  of  a  few  per  cent  in  the  number  of  free  electrons.  However,  mul¬ 
tiple  ionization,  including  Auger  effect,  was  already  included  in  the 
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experiments  of  Schram  et  ci,  ^ince  the/  measured  electric  current  rather 
than  counting  ionization  events. 

To  estimate  the  time  to  form  the  residual  ionization,  we  must  know  some¬ 
thing  about,  the  spectrum  of  ener/v  «  of  the  secondary  electrons  (i.e., 
primary  ionization) .  There  appears  to  be  no  experimental  data  on  this  spec¬ 
trum.  Theoretically,  it  is  expected  to  be  approximately  of  the  form  dw/w*, 
as  in  Rutherford  scattering.  Atomic  binding  will  provide  a  cutoff  of  the 
divergence  at  w+0.  We  shall  therefore  assume  that  the  number  dn  of  secon¬ 
dary  electrons  having  energy  in  dw  is: 


dn  = 


A  dw 
w„2+  w2 


To  normalize  this  form  to  unity,  we  integrate  from  w^O  to  w=ff/2,  where  W 
is  the  energy  of  the  primary  electron;  we  define  the  secondary  electron  to 
be  that  one  of  the  two  outgoing  electrons  which  has  the  smaller  energy: 

W/2 


tan  = 


arctan 


/  K  \  a  A_ 

V2v  2  *= 


In  the  iast  form  here,  we  have  taken  advantage  of  the  fact  that  W/2u,  is 
(usually)  a  very  large  number. 

From  Equations  (76)and  (77),  we  can  calculate  the  mean  energy  of  thr 

secondary  electrons: 

W/2 


«  *  £ if  Zp-r  - 

‘  j  5V*1  ”  \  4“t'  / 


This  result  can  be  used  to  determine  an  approximate  value  for  the  constant 
wa.  We  hax'e  seen  above  that  for  each  free  electron  produced,  the  primary 
electron  loses,  on  the  average,  86  ev.  Of  this,  we  estimate  20  ev  is 
expended  for  the  ionization  potential  3nd  another  10  ev  left  in  excitation 
of  atoms  (or  molecules) .  Thus  the  mean  kinetic  energy  w  of  the  secondary 
electrons  should  be  about  56  ev.  Taking  W  =  10e  ev  in  Equation  (78), 
one  then  finds: 


We  see  that  equal  amounts  of  secondary  energy  are  contained  in 
each  decade.  Since  the  amount  of  secondary  ionization  is  approximately 
proportional  to  the  kinetic  energy  available,  the  secondary  ionization 
also  will  come  equally  from  each  decade. 

Figure  (4)  is  a  graph  of  o^v  for  electrons  on  air  atoms.  In 
plotting  this  graph  we  have  used  experimental  data  directly,  rather  than 
Eq.  (71),  which  is  not  very  accurate  for  energies  less  than  about  103  ev. 


The  initial  rate  R0  of  production  of  ionization  by  secondary  elec¬ 
trons,  per  secondary  electron,  is 


R0  =  N 


f  o.v  dn  =  N  f  — ■  -a d(£n  w) 
J  i  J  ^  w02+w2 


A  numerical  evaluation  of  this  integral  leads  to  the  result 


R0  =  (1.04  x  10'8  cm3/sec. )  N  . 


This  rate  of  production  of  ionization  will  not  be  maintained  very 
long.  Fig.  (5)  is  a  graph  of  the  integrand  in  Eq.  (82),  i.e.  the  factors 
multiplying  d(£n  w) .  One  sees  that  most  of  the  integral  (5/6  of  it,  in 
fact)  comes  from  electrons  in  the  energy  decades  below  about  250  ev. 
According  to  Eq.  (SI),  or  the  table  following  it,  only  about  50%  of  the 
total  secondary  ionization  arises  from  secondary  electrons  in  this  energy 
range.  The  initial  ionization  rate  therefore  cannot  persist  for  a  time 
longer  than 


0.3  x  1.52 
|  1.04xi0~®N 


5.3  :c  U)7 


In  fact,  the  rate  should  have  fallen  by  a  factor  4  or  5  by  this  time, 
since  the  secondaries  with  energy  less  than  250  ev  will  have  been 
exhausted. 


Figure  4.  Graph  of  o{v  as  a  function  of  w  for  electrons  on  air  atoms. 
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I , (T)  =  0.160  £n 


(i.  — g.  -r)  . 

\  1.54x10'/ 


(95) 


Thus  I2(T)  increases  only  slowly  with  time,  and  reaches  half  its  final 
value  at 


NT.  -  10"2  NT 
*/2  max 


2xl09  sec/cm3  . 


(96) 


In  EMP,  one  is  often  interested  in  the  amount  of  ionization  pro¬ 
duced  in  a  time  of  10“8  second  after  loss  of  a  given  amount  of  energy  by 
the  primary  (Compton  recoil  electron).  Per  34  ev  lost  by  the  primary, 
the  number  of  ion  pairs  formed  by  this  time  is 


N. 

IP 


34 

86 


[,  .  0.160  f.„(l  .  -r^nr)]  • 


(9?) 


N. 

i? 


This  number  is  given  in  Table  5  as  a  function  of  altitude.  One  sees 
that  the  effect  of  the  time  lag  in  formation  of  the  secondary  ionization 
is  an  important  effect  at  altitudes  above  40  km  or  so. 


♦■itude,  km  = 


Table  5 

0  10  20  30  40  50 

1.00  0.98  0.90  0.79  0.69  0.6) 


Electron  Equilibration 
Once  the  secondary,  tertiary,  etc.,  electrons  have  fallen  below 
about  15  ev,  they  can  no  longer  produce  additional  ionization.  However, 
they  still  lose  energy  in  collisions  with  air  molecules.  If  an  electric 
field  0  is  present,  the  electrons  gain  energy  from  it.  After  some  time, 
the  electrons  reach  an  equilibrium  distribution  in  energy,  in  which  energy 
lost  to  air  molecules  is  balanced  by  energy  gained  from  the  electric  field. 
Tills  equilibrium  distribution  depends  only  on  the  ratio  E/p,  where  p  is  the 
air  pressure,  if  the  air  temperature  is  assumed  constant. 
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In  EhIP  calculations,  it  is  usually  assumed  that  the  electrons 
acquire  the  equilibrium  distribution  immediately  after  birth,  and  also 
follow  changes  in  E  (or  E/p)  without  time  lag.  Under  this  assumption, 
such  important  quantities  as  the  electron  attachment  rate  (to  02),  the 
electron  drift  velocity  or  mobility,  and  the  cascading  rate  or  Townsend 
coefficient,  depend  only  on  the  instantaneous  value  of  E/p.  Conveniently, 
experimental  data  for  these  quantities  is  usually  obtained  as  functions 
of  E/p.  Thus,  in  this  approximation,  it  is  not  necessary  to  know  the 
details  of  the  electron  energy  distribution. 

If  we  now  wish  to  take  into  account  the  deviations  from  equili¬ 
brium,  including  time  lags,  we  shall  have  to  introduce  parameters  which 
characterize  the  electron  energy  distribution.  The  simplest  such  charac¬ 
terization  is  to  use  only  the  mean  energy  of  the  distribution.  Actually, 
the  parameter  Ug  of  electron  swarm  theory  is  2/3  of  the  mean  energy, 
in  analogy  to  temperature,  although  the  energy  distribution  is  not  Max¬ 
wellian  if  an  electric  field  is  present. 

We  have  seen  in  Sec.  S  that  about  70%  of  the  electrons  have  energy- 
less  than  16  ev  at  birth.  Most  of  the  other  30%  drop  rapidly  (in  a 
time  given  by  Eq.  (84)  belGw  16  ev  by  producing  further  ionization.  There 
fore,  it  seems  reasonable  to  start  the  electrons  off,  in  equilibration 
calculations,  with  a  mean  energy  of  about  8  ev,  or  Ug  about  5  ev. 

Baum  (Ref.  7)  has  studied  the  relaxation  of  U  ,  using  data  on 
momentum  and  energy  transfer  collision  frequencies  provided  by  A. V. Phelps 
We  may  use  Baum's  results  to  appraise  the  importance  of  the  finite  relax¬ 
ation  time  for  BMP  calculations.  We  take  as  a  typical  electric  field 
I  esu  =  3xlGu  volt/meter.  Table  6  then  shows,  for  various  altitudes, 
the  equilibrium  value  of  and  the  time  to  reach  equilibrium,  starting 

from  U  =  5  ev.  One  sees  that,  for  the  electric  field  assumed,  the  times 
are  all  less  than  10~9  second.  It  therefore  appears  that  the  finite 
relaxation  time  cannot  have  a  drastic  effect  (by  altering  the  effective 


electron  mobility)  on  the  peak  electric  field,  unless  the  computed  elec¬ 
tric  field  changes  appreciably  in  10  5  seconds. 

Table  6 


altitude,  km 

=  0 

10 

20 

30 

40 

ue,  ev 

=  0.17 

0.35 

1.0 

1.6 

time,  sec 

=  7x10"10 

9x10" 10 

3x10" 50 

U* 

X 

— i 

o 

1 

o 

M0" 

Note  that  at  40  km,  the  equilibrium  Ue  is  high  enough  that,  for  a 
Maxwell  distribution  of  electron  energies,  cascading  would  occur  rapidly. 
However,  for  high  Ue,  the  deviations  from  the  Maxwellian  distribution 
are  large.  To  compute  the  cascading  rate  correctly,  one  would  need  more 
information  (than  Ue)  concerning  the  electron  distribution.  However,  there 
does  not  appear  to  be  enough  basic  cross-section  data  available  (currently) 
to  permit  an  adequate  treatment  of  the  detailed  distribution.  Thus  the 
cascading  rate  can  be  computed  only  for  the  equilibrium  case,  where  it  is 
known  from  the  Townsend  coefficient,  which  has  been  determined  experi¬ 
mentally  as  a  function  of  E/p.  The  determination  of  the  cascading  rate 
from  the  data  on  the  Townsend  coefficient  is  discussed  in  Ref.  8. 

At  low  altitude,  where  cascading  normally  is  not  important,  one 
could  devise  a  non-equilibrium  correction  by  using  Baum's  relaxation  times. 
It  would  be  important  to  do  this  oniy  if  rise  times  of  the  electric  field 
are  as  short  as  10-?  second. 


33 


References 

(1)  See  the  review  article  by  Bethe  and  Ashkin  in  "Experimental 
Nuclear  Physics",  E.  Segre,  editor,  Vol.I,  fe'iley,  1953;  p.253 

(2)  See  the  discussion  in  Mott  and  Massey,  "The  Theory  of  Atomic 
Collisions",  Third  Edition,  Oxford,  1965;  p.467 

(3)  See  Longmire,  "Elementary  Plasma  Physics",  Wiley-Interscience, 
1965;  p. 199 

(4)  Marshall  and  Ward,  Can.  J.  Research,  A15,  39  (1937) 

(5)  H.A.  Sethe,  Ann  Physik,  £,  325  (1930);  also  2.  Physik,  76, 

293  (1938) 

(6)  B.L.  Schraro,  F.J.  Dellecr,  M.J.  Van  der  Kiel  and  J.Kistenacher, 
Physica,  31^,  94  (1965) 

(7)  ARYL  EMP  2-1,  Electromagnetic  Pulse  Theoretical  Notes, 

Volume  1,  Note  12,  by  C.E.Saum. 

(8)  MRC-R-  3  ,  ''Development  of  the  CHAP  EMP  Code  (U)"  by  H.J. 
longley  and  C.L. Longmire,  in  publication.  Unclassified 


j  w-.-  ~  CZ 


